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Abstract 


The  spatial  frequency  response  of  an  imaging  system,  known  as  the  Modulation  Transfer 
Function  (MTF),  is  a  primary  image  quality  metric  that  is  commonly  measured  with  a 
sine  wave  target.  The  FBI  certification  program  for  commercial  fingerprint  capture 
devices,  which  MITRE  actively  supports,  has  an  MTF  requirement.  In  some  cases, 
however,  a  square  wave  ("bar  target")  must  be  used  in  testing,  which  results  in  a  similar 
quantity  called  the  Contrast  Transfer  Function  (CTF).  This  document  reports  on  an 
investigation  of  the  mathematical  relationship  between  the  MTF  and  CTF,  methods  for 
converting  between  the  two,  and  derives  an  equivalent  CTF  from  the  given  spec  MTF,  for 
use  in  the  FBI  certification  program.  The  methodology  presented  is  applicable  to  the 
general  case,  i.e.,  whenever  conversion  between  the  MTF  and  CTF  of  an  imaging  system 
is  needed. 
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Section  1 


Introduction 

1.1  Background 

The  Federal  Bureau  of  Investigation  (FBI)  recently  fielded  an  Integrated  Automated  Fingerprint 
Identification  System  (lAFIS)  that  greatly  enhances  the  FBI’s  identification  services  to  other 
federal  agencies,  and  to  local,  state,  and  regional  law  enforcement  agencies.  lAFIS  is  a  large, 
complex  information  system  that  utilizes  off-the-shelf  and  advanced  technologies.  Agencies 
interacting  with  lAFIS  require  some  components  of  these  technologies  to  maintain 
interoperability. 

Among  other  components,  lAFIS  incorporates  fingerprint  scanners  to  convert  inked  fingerprint 
cards,  hardcopy  latent  fingerprints,  and  live  subjects'  finger  scans  (livescan)  into  digital  data,  at 
either  500  pixels  per  inch  (ppi)  or  1000  ppi,  with  output  at  8  bits  per  pixel  (256  gray  levels). 
MITRE  is  supporting  the  FBI's  certification  process  for  vendor  fingerprint  products,  whose 
foundation  is  the  lAFIS  Image  Quality  Specifications  (IQS)  (FBI,  1999).  The  IQS  defines  the 
quantitative  image  quality  requirements  for  fingerprint  scanners  and  printers.  Products  are 
certified  by  a  process  whereby  the  vendor  generates  test  data  by  following  a  set  of  established 
test  procedures  (FBI,  1995)  and  MITRE  analyzes  the  test  data  for  conformance  and  compliance 
with  the  specification  requirements.  This  successful  certification  process  has  helped  to  ensure  a 
standard  level  of  image  quality  for  fingerprint  comparisons,  classification,  feature  detection,  and 
automated  search  reliability,  and  supports  systems  interoperability  between  national,  state,  and 
local  law  enforcement  agencies. 

One  of  the  critical  components  of  the  IQS  test  suite  is  the  testing  performed  to  ascertain  a 
fingerprint  scanner's  spatial  frequency  response,  known  as  the  Modulation  Transfer  Eunction 
(MTE).  Eundamentally,  the  MTE  is  a  measure  of  the  contrast  transmission  capabilities  of  an 
imaging  system,  as  a  function  of  spatial  frequency.  The  MTE  has  proven  to  be  an  excellent 
measure  of  the  sharpness  and  detail  rendition  of  an  imaging  system  (Tannas,  1985).  Eor 
fingerprint  scanners,  it  indicates  the  devices’  capability  to  capture  fingerprints  at  a  quality  level 
that  is  needed  for  successful  latent  image  comparison,  as  well  as  for  inked  or  livescan  fingerprint 
search/matching  and  examination.  MTE  assessment  has  successfully  quantified  the  elusive  term 
‘image  quality’  across  fingerprint  capture  devices  of  many  different  designs. 

A  sine  wave  target  is  specified  for  use  in  determining  the  MTE  in  IQS  testing.  A  sine  wave  is  the 
fundamental  pattern  for  directly  determining  the  spatial  frequency  response  of  a  linear  system;  a 
pure  sine  wave  input  to  a  linear,  spatially  invariant  system  is  output  as  a  pure  sine  wave  of  the 
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same  frequency,  generally  with  reduced  amplitude\  It  has  a  high  signal-to-noise  ratio,  and  the 
fact  that  each  spatial  frequency  is  individually  imaged  has  been  found  to  be  very  useful  in  fault 
detection/isolation  of  scanner  problems.  In  addition,  the  maximum  value  MTF  is  the  IQS- 
specified  quantity  of  interest  and  a  sine  wave  target  readily  lends  itself  to  such  a  measurement, 
because  the  single  sine  period  in  each  frequency  pattern  that  exhibits  the  highest  computed 
modulation  can  be  readily  identified.  Finally,  the  target  is  available  as  a  commercial  off-the-shelf 
product^;  a  typical  sine  wave  target  is  shown  in  Figure  1 . 


Figure  1.  Typical  Sine  Wave  Target 


'  Other  targets  can  be  used  to  measure  the  MTF,  e.g.,  a  narrow  line  target,  point  target,  random  pattern,  or  sharp 
edge.  In  these  cases  the  target  is  broken  up  into  its  sine  wave  components,  via  Fourier  analysis,  to  compute  the  sine 
wave  MTF.  Each  target  type  has  its  own  set  of  advantages  and  disadvantages. 

^  Three  vendors  (first  is  manufacturer);  http://www.sinepatterns.com,  http://www.edmundoptics.com, 
http://www.jmloptical.com 
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For  IQS  testing,  the  sine  wave  target  is  scanned  and  the  resulting  digital  image  is  processed 
through  a  MITRE-developed  computer  algorithm  which  computes  the  scanner’s  MTF  (MITRE, 
2001).  The  essence  of  the  technique  is  illustrated  in  Eigure  2  for  a  single  sine  wave  frequency:  the 
gray  levels  corresponding  to  the  peak  and  valley  in  the  sine  wave  image  are  identified,  these  gray 
levels  are  converted  to  target  space  reflectance  values  by  traversing  the  scanner's  input/output 
response  curve,  the  image  modulation  is  computed,  and  divided  by  the  target  modulation  to 
obtain  a  single  point  on  the  MTE  curve.  The  process  is  then  repeated  for  a  series  of  sine  wave 
frequencies  until  the  complete  MTE  curve  is  generated. 

A  complication  in  applying  the  MTE  concept  to  image  scanners  results  from  the  fact  that 
scanners  perform  discrete  sampling  over  finite  areas.  A  scanner  samples  a  continuous,  analog 
image  input  at  discrete  locations  over  discrete  time  intervals,  integrating  all  incident  light  energy 
over  a  given  time  interval  and  within  a  small  area  around  each  sampling  point,  where  the  small 
area  is  the  individual  detector  element  ("pixel")  area.  One  of  the  consequences  of  this  discrete 
sampling  by  finite- sized  detector  elements  is  that  the  scanner  acts  as  a  space  variant  system, 
which  implies  that  the  scanner  MTE  can  be  multivalued  (Park,  et.al,  1984;  Eeltz,  et.al.,  1990).  In 
the  context  of  the  sine  wave  target,  this  space  variance  means  that  the  phase  between  this  target 
and  the  scanner's  detector  array  plays  a  role  in  the  actual  MTE  achieved.  Here,  phase  refers  to 
the  location  of  the  detector  pixel  within  a  sine  wave  period,  at  the  time  the  pixel  is  collecting  light 
energy.  As  the  target-detector  phase  changes,  the  scanner  MTE  can  change,  within  limits.  In 
practice,  one  measures  the  minimum,  maximum,  or  average  MTE.  It  has  been  established 
procedure  in  IQS  testing  to  use  the  maximum  measured  MTE  (after  noise  filtering)  for  comparing 
to  the  spec  MTE  for  compliance  verification  purposes. 
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Figure  2.  Path  of  Single  Sine  Wave  to  Generate  MTF  Data  Point 
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1.2  Problem  Statement 


There  are  two  fundamentally  different  types  of  fingerprint  scanners  certified  for  use  in  the  FBI 
lAFIS  program.  In  one  type,  inked  fingerprints  on  card  stock  ("ten-print  card")  or  latent  prints 
on  paper  are  scanned.  Typically,  this  scanner  is  a  general  purpose,  Commercial-Off-The-Shelf 
(COTS)  or  modified  COTS  flatbed,  desktop  scanner  operating  in  reflection  mode,  with  a  linear 
Charge-Coupled  Device  (CCD)  scanning  array  detector  and  associated  signal  processing  to 
produce  the  digital  output  image.  The  common  imaging  mode  utilized  in  these  scanners  is  always 
compatible  with  scanning  a  reflective  sine  wave  target. 

The  other  type  of  scanner,  here  called  a  livescanner.  directly  scans  a  person's  finger  and  is 
designed  solely  for  that  purpose.  In  the  most  common  livescanner  design  a  subject's  finger  is 
pressed  against  ("plain  impression")  or  rolled  across  ("rolled  print")  the  surface  of  a  glass  prism. 
The  image  is  created  by  the  difference  in  index  of  refraction  (n)  between  skin  and  air,  utilizing  the 
optical  imaging  principle  of  "total  internal  reflection"  (TIR),  also  called  "frustrated  internal 
reflection"  (FIR)  (Smith,  1966).  Within  certain  defined  angles  of  incident  illumination,  light 
traveling  from  the  prism  surface  to  finger  ridges  in  contact  with  the  glass  will  be  absorbed  by  the 
skin  (ngiass  =  nski„),  while  light  traveling  from  the  prism  surface  to  finger  valleys,  where  a  glass-air 
gap  exists  (ngiass  >  nair),  will  be  reflected  back  into  the  glass.  The  reflected  light  is  detected  by  a 
linear  or  two-dimensional  CCD  array  to  produce  the  digital  output  image. 

FIR  imaging  in  livescanners  depends  on  the  three-dimensionality  of  skin,  which  is  either  in 
contact  or  not  in  contact  with  a  glass  surface,  and  the  relative  indices  of  refraction  between  glass, 
air,  and  skin.  Many  livescanners  cannot  adequately  image  a  two-dimensional,  continuous  gray 
tone  sine  wave  target;  the  output  image  suffers  from  low  signal-to-noise  ratio,  low  contrast,  and 
very  few  gray  levels,  making  it  impractical  to  compute  a  valid  MTF.  It  has  been  found  through 
experimentation,  however,  that  livescanners  can  adequately  image  reflective  or  transmissive  bar 
targets,  mainly  because  a  bar  target  is  extremely  high  contrast  with  only  two  gray  levels,  black 
and  white,  and  it  can  be  easily  fabricated  on  a  glass  substrate. 

The  fundamental  problem  addressed  in  this  document  can  be  stated  as  follows: 

The  FBI  image  quality  specification  specifies  the  minimally  acceptable  values  for  a  scanner  MTF 
when  it  is  measured  using  a  sine  wave  target.  Since  a  bar  target  must  be  used  in  some  product 
certification  cases,  the  problem  is  to  answer  the  questions:  (1)  is  there  a  relationship  between  the 
spatial  frequency  response  derived  from  a  bar  target  and  the  MTF  derived  from  a  sine  wave 
target?,  (2)  what  is  this  relationship?,  and  (3)  what  is  the  best  approach  for  handling  a  bar  target 
derived  spatial  frequency  response,  given  that  the  specification  is  for  a  sine  wave  target  derived 
MTF? 
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Section  2 


Problem  Solution  -  Synopsis  of  Results 

The  analyses  and  theory  presented  over  the  years  in  technical  literature  confirms  that,  for  the 
same  imaging  system,  the  bar  target  derived  spatial  frequency  response  and  sine  target  derived 
MTF  are  not  equal  to  each  other,  but,  as  will  be  seen,  there  is  a  definite  mathematical  relationship 
between  them.  The  difference  between  the  two  quantities  is  often  highlighted  by  denoting  the 
spatial  frequency  response  obtained  from  a  bar  target  as  the  Contrast  Transfer  Function  (CTF)  or 
square  wave  response,  to  distinguish  it  from  the  term  'MTF',  which  is  reserved  for  the  sine  wave 
transfer  function.  This  terminology  distinction  between  bar  target  CTF  and  sine  target  MTF  is 
used  in  the  remainder  of  this  document. 

In  solving  the  problem  at  hand,  where  a  sine  wave  target  MTF  is  specified  but  the  imaging 
system  under  test  may  require  use  of  a  bar  target,  one  can  either: 

(1)  convert  the  spec  sine  wave  MTF  to  the  equivalent  bar  CTF,  which  then 
becomes  the  spec  CTF  when  bar  targets  are  used, 

or, 

(2)  convert  each  CTF  measured  with  a  bar  target  to  its  equivalent  sine  wave  MTF,  for 
comparison  to  the  original  spec  sine  wave  MTF. 

The  first  approach  is  more  straightforward  to  apply  in  the  testing  environment,  with  less  chance 
of  computational  error,  because  the  modulations  of  a  sample  are  directly  computed  from  the  bar 
target  image  and  compared  to  the  'spec'  bar  CTF,  with  no  further  manipulation  or  conversion  of 
sample  data.  Therefore,  the  first  approach  is  the  solution  approach  utilized. 

The  method  of  construction  of  the  spec  CTFs  are  discussed  in  detail  in  Section  3;  here  we 
present  the  final  results  in  terms  of  the  proposed  spec  CTFs  for  IQS  certification  compliance 
verification.  Table  1  gives  the  proposed  spec  CTFs  for  a  500  ppi  scanner  and  for  a  1000  ppi 
scanner,  which  are  equivalent  to  the  corresponding  IQS  Appendix  F  spec  sine  wave  MTFs,  also 
given  in  Table  1  (see  also  Figure  6).  The  following  points  should  be  noted  with  respect  to 
Table  1: 

•  The  CTF  values  in  Table  1  are  unnormalized  and  assume  a  target  modulation  of  1.0  at  all 
frequencies. 
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•  In  IQS  testing,  the  bar  target  frequencies  cannot  be  predefined  because  the  exact  target,  either 
COTS  or  special  design,  is  a  vendor  choice.  The  actual  CTF  is  moderately  dependent  on  the 
number  of  bars  at  a  given  frequency;  Table  1  values  are  the  average  of  5 -bar,  15 -bar,  and  50-bar 
CTFs.  A  bar  target  with  less  than  5  bars  per  frequency  should  never  be  used  for  CTF 
measurement^. 

•  For  a  bar  target  containing  frequencies  not  given  in  Table  1,  the  following  equations  can  be  used 
to  obtain  the  CTF  modulation  values  at  those  other  frequencies. 

500  ppi,  f  =  1.0  to  10.0  cy/mm, 

CTF=3.04105E-04*f"-7.99095E-02*f  -1-1.02774  (1) 

1000  ppi,  f  =  1.0  to  20.0  cy/mm, 

CTE  =  -1.85487E -  05 *  f^  +  1.41666E  -  03  *  f"  -  5.73701E  -  02  *  f  +  1.01341  (2) 


3 

A  3-bar  target,  such  as  the  USAF  1951  tribar  target,  is  not  appropriate  for  CTF  measurement  because  the  peak  and 
valley  image  amplitudes  are  erratic,  only  the  center  bar  and  two  adjacent  spaces  are  available  for  modulation 
calculation,  and  because  obtaining  good  phasing  between  the  target  and  scanner  detector  pixels  is  difficult.  A 
minimum  of  a  5-bar  target  is  required  for  stability,  such  as  the  NBS  lOlOA  target  (a.k.a.,  NBS  1963A  target, 
ANSl/lSO  Test  Chart  #2);  better  is  a  15 -bar  target,  such  as  the  T90  target. 
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Table  1.  IQS  Spec  MTFs  and  Corresponding,  Proposed  Spec  CTFs 


Frequency 

(cy/mm) 


0.5 


1.0 


1.5 


2.0 


2.5 


3.0 


3.5 


4.0 


.5 


5.0 


5.5 


6.0 


6.5 


7.0 


7.5 


8.0 


8.5 


9.0 


9.5 


10.0 


10.5 


11.0 


11.5 


12.0 


12.5 


13.0 


13.5 


14.0 


14.5 


15.0 


15.5 


16.0 


16.5 


170 


17.5 


18.0 


18.5 


19.0 


19.5 


20.0 


MTF 
500  ppi 


CTF 
500  ppi 


MTF 
1000  ppi 


0.210 


CTF 
1000  ppi 


0.979 


0.957 


0.930 


0.904 


0.879 


0.854 


0.829 


0.805 


0.782 


0.760 


0.738 


0.716 


0.695 


0.675 


0.655 


0.636 


0.617 


0.598 


0.580 


0.563 


0.546 


0.529 


0.513 


0.497 


0.481 


0.466 


0.451 


0.437 


0.423 


0.409 


0.395 


0.382 


0.369 


0.356 


0.344 


0.332 


0.319 


0.308 


0.296 


0.284 
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Section  3 


Technical  Details  of  Solution 

A  straightforward  series  expansion  formula  derived  nearly  half  a  century  ago  (Coltman,  1954) 
will  quite  accurately  convert  a  linear  imaging  system's  sine  wave  target  MTF  to  its  equivalent  bar 
target  CTF,  or  vice  versa,  under  certain  conditions.  These  conditions  include  the  assumption  that 
the  imaging  system  is  analog,  e.g.,  a  photographic  film  camera,  and  the  assumption  that  the 
number  of  target  bars  is  infinite,  i.e.,  the  bar  target  is  a  true  square  wave  of  infinite  extent.  The 
Coltman  formula  and  its  application  is  described  in  more  detail  in  Section  3.2,  but  note  that  the 
implicit  assumptions  in  this  formula  do  not  hold  for  the  case  of  interest  here.  Specifically,  a 
digital  scanner  is  a  discrete  sampling  system,  not  a  continuous  analog  system,  and  a  real  bar  target 
used  to  measure  the  CTF  contains  a  finite  number  of  bars,  sometimes  as  few  as  5  bars  per 
frequency.  Under  these  conditions,  and  using  a  curve  extrapolation  strategy,  the  Coltman 
formula  can  still  give  a  reasonably  accurate  CTF,  as  discussed  in  Section  3.2. 

However,  we  seek  a  method  for  converting  the  MTF  to  the  CTF  which  directly  takes  into 
account  the  discrete  sampling  of  a  digital  scanner,  the  finite  number  of  bars  in  a  bar  target,  and 
results  in  a  high  accuracy  CTF,  which  is  a  true-equivalent  of  the  FBI's  IQS  specified  sine  wave 
MTF.  The  method  developed  in  Section  3.1  meets  these  conditions. 


3.1  Solution  Approach 

3.1.1  Methodology 

The  following  procedure  was  developed  for  converting  a  discrete  sampling  imaging  system's 
MTF  to  its  equivalent  CTF  for  a  given  number  of  target  bars: 

(a)  the  given  system  MTF  is  divided  by  the  Fourier  transform  of  the  sampling  pixel, 

(b)  the  result  of  (a)  is  multiplied  by  the  Fourier  transform  of  the  bar  target, 

(c)  the  result  of  (b)  is  Fourier  transformed, 

(d)  the  result  of  (c)  is  convolved  with  the  sampling  pixel,  and 

(e)  the  maximum  and  minimum  values  from  the  result  of  (d)  are  located,  from  which  the  CTF 
modulation  values  are  computed. 

Steps  a-e  are  illustrated  in  the  following  math-graphic,  where  F  denotes  Fourier  transformation 
and  0  denotes  convolution.  Steps  a-e  are  individually  discussed  in  detail  in  the  following. 
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bar  tgt  pixel 


eval 

max  &  min 


Step  (a)  -  For  an  analog  imaging  system,  the  MTF  ean  be  measured  up  to  the  frequeney  at  whieh 
the  MTF  goes  to  zero,  which  is  commonly  called  the  'cut-off  spatial  frequency  in  the  optical 
field,  equivalent  to  the  'band-limit'  frequency  in  electrical  engineering  parlance.  This  step  is 
performed  in  order  to  obtain  the  imaging  system  MTF  out  to  its  cut-off  frequency,  i.e.,  in  analog 
space,  before  the  image  is  sampled  by  the  discrete  detector  pixels.  For  the  case  at  hand,  the  given 
system  MTF  is  the  FBI  IQS  spec  MTF,  which  is  truncated  at  the  sampling  pixel's  Nyquist 
frequency,  where  the  modulation  is  non-zero'^.  In  equational  form,  the  procedure  is  given  by: 


M(f)  = 


M.(f) 

sinc(7tfw) 


(3) 


where, 

M(f)  =  imaging  system  MTF  before  discrete  detector  sampling 
Ms(f)=  total  imaging  system  MTF  =  IQS  spec  MTF 
sinc(7Cfw)  =  sin(7tfw)  /  Ttfw 
f  =  spatial  frequency  in  cy/mm 

w  =  detector  pixel  width:  0.0508  mm  for  500  ppi,  0.0254  mm  for  1000  ppi 

Since  the  spec  MTF  is  only  defined  up  to  Nyquist  frequency,  extrapolation  must  be  used  to 
extend  M(f)  until  it  reaches  zero  modulation.  A  reasonable  approach  is  to  monotonically  and 
smoothly  decrease  modulation  beyond  Nyquist  until  it  reaches  zero.  This  was  obtained  by  fixing 
the  zero  modulation  at  four  times  the  Nyquist  frequency  and  fitting  a  smooth  curve  between 
Nyquist  and  four  times  Nyquist  frequency.  This  approach  makes  sense  because  M(f)  is 
essentially  the  scanner's  optics  MTF  and  an  optics  MTF  is  everywhere  positive  and  smoothly. 


The  Nyquist  frequency  equals  the  reciprocal  of  twice  the  pixel  pitch,  where  pixel  pitch  is  the  pixel-center-to- 
adjacent-pixel-center  distance.  In  a  real  scanner,  the  pixel  pitch  may  be  greater  than  or  equal  to  the  pixel  width.  For 
example,  for  500  ppi  sampling  output  with  a  scanner  having  a  true  optical  resolution  of  1000  ppi,  the  pixel  pitch  = 
1/500  =  0.02"  while  the  pixel  width  is  0.01".  For  a  discrete  sampling  system,  valid  MTF  measurements  cannot  be 
obtained  above  Nyquist  frequency.  The  spec  MTF  has  a  value  of  0.2  modulation  at  Nyquist  frequency  but  a  scanner 
could  have  a  value  as  high  as  -0.7  at  Nyquist.  In  either  case,  the  MTF  curve  is  truncated  at  Nyquist,  making  it 
computationally  intractible. 
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monotonically  decreases  to  zero  frequency^.  Experiments  run  demonstrate  that  CTF  results  up 
to  Nyquist  frequency  are  not  greatly  affected  by  any  reasonable  choice  of  M(f)  values  above 
Nyquist  frequency. 

For  computational  efficiency,  a  curve  fit  to  the  entire  M(f)  was  performed,  with  the  following 
results  for  500  ppi  and  1000  ppi  MTFs;  also  plotted  in  Figure  3. 

500  ppi,  f  <  10.0  cy/mm: 

M(f)=4.8597E-07*f  -  1.856125E05  *f®  +  2.846349E- 04*f^  -  2.300554E03*F‘ 

(4) 

+  1.060994E2*f^-  2.276723E  2*f^  -  7.703449E  2  *f  +  1.0000090 
500  ppi,  f  >  10.0  cy/mm: 

M(f)=exp(-9.991357&02*f'“)  (5) 

1000  ppi,  f  <  20.0  cy/mm: 

M(f)  =  -2.99661568E8*f^  +  3.04795483E6 *f'‘  -  1.35961019E4*f^ 

+  3.9614024E3*f^  -  7.76516183E2*f  +  1.00001941  (6) 

1000  ppi,  20.0  <  f  <  80.0  cy/mm: 

M(f)  =  -9.90555556E7*f^  +  2.1147619E4*f^-  .O184142778E0=f  +  0.627922857  (7) 


Steps  (b),  (c)  -  These  steps  are  conveniently  performed  by  utilizing  a  closed  form  equation  for 
computing  the  image  profile  of  an  N-bar  target,  given  a  system  MTF,  that  has  been  reported  in 
the  literature  (Barakat  and  Ferman,  1967,  Equation  (10)).  With  a  change  of  notation,  Barakat's 
Equation  (10)  becomes: 


cos(27tx/)sin(47tAT/)M(/) 

/cos(27iL/) 


(8) 


where, 

j(x)  =  image  intensity  profile  as  a  function  of  distance  x,  prior  to  pixel  sampling 
N  =  number  of  bars  in  bar  target,  at  specific  frequency 
F  =  half-bar  width 

fc  =  cut-off  frequency  of  M(f):  40  cy/mm  for  500  ppi,  80  cy/mm  for  1000  ppi 


^  For  the  special  cases,  not  considered  here,  of  optics  with  large  aberrations,  large  defocus,  or  a  central  obscuration 
(apodized  aperture),  the  optics  MTF  still  goes  to  zero,  but  not  necessarily  monotonically.  Also  note  that,  at  first 
glance,  it  might  seem  that  the  spec  MTF  could  be  extrapolated  first,  and  then  divided  by  sinc(.);  but  this  approach 
would  not  give  a  smooth  result  because  sinc(.)  oscillates  and  goes  through  zero  modulation,  at  frequencies  above 
2xNyquist. 
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Step  (d)  -  Equation  (8)  is  convolved  with  a  rectangle  function  representing  the  pixel  sampling 
performed  in  a  scanner,  to  obtain  the  final  output  image  profile: 


i(x)  =  r(t)  j(x  -  t)dt 


(9) 


where, 

r(t)  =  rectangle  function,  width  =  0.0508  mm  for  500  ppi, 
width  =  0.254  mm  for  1000  ppi 

Step  (e)  -  Finally,  the  maximum  and  minimum  values  of  the  total  imaging  system's  bar  image 
profile  are  identified  from  i(x).  This  is  equivalent  to  the  integral  of  r(x)j(x)  with  r(x)  centered  at  a 
valley  in  j(x)  and  then  centered  at  a  peak  in  j(x).  An  alternative  method  is  to  constrain  all  r(x) 
locations  to  be  integer  multiples  of  the  pixel  width,  and  then  determine  which  two  r(x)  locations 
result  in  maximum  modulation.  In  this  case  the  maximum  and  minimum  may  not  be  centered  on 
the  peak  and/or  valley  of  j(x),  respectively.  In  modeling  an  actual  scanner's  performance  this 
second  method  may  be  more  realistic,  but  for  identifying  the  CTF  that  most  closely  corresponds 
to  a  spec  MTF,  the  first  method  makes  more  sense. 


Spatial  Frequency  (cy/mm) 

Figure  3.  Curve  Fits  to:  M(f)  =  Spec  MTF  /  Sinc(7tfw) 
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3.1.2  Computer  Program 


A  computer  program  was  developed  in  Fortran-90  to  perform  steps  b-e,  whieh  ineorporates  the 
separately  eomputed  eurve  fits  given  in  Equations  (4-7).  Some  features  of  this  eomputer 

program  are  as  follows: 

•  Caleulations  are  performed  in  double  preeision  mode. 

•  Integrations  are  performed  via  Simpson's  rule  with  a  large  number  of  small  integration 
inerements  to  aeeount  for  the  high  frequeney  oseillations  of  the  trig  funetions. 

•  Evaluation  is  only  performed  over  1.25  periods  of  the  bar  target,  eentered  on  the  middle  bar, 
beeause  the  highest  peak  and  lowest  valley  always  oeeur  at  and  adjaeent  to  the  middle  bar  in  this 
symmetrieal  model. 

•  The  input  bar  target  is  assumed  to  have  a  modulation  of  1.0. 

A  number  of  eode  verifieation  tests  were  run;  the  major  tests  were  as  follows: 

•  Identified  required  integration  inerements  by  running  series  of  inerement  sizes  and  observing 
where  stable  output  oeeurred. 

•  Proeessed  a  50-eyele  sine  wave  target  with  0.0508  mm  and  0.0254  mm  sampling  pixel  widths 
and  with  the  500  ppi  and  1000  ppi  spec  MTEs,  respectively.  The  output  MTE  should  equal  the 
original  spec  MTE  in  each  case,  whieh  they  did  to  within  less  than  0.4%  for  500  ppi  ease  and  to 
within  less  than  0.7%  for  the  1000  ppi  ease;  these  results  are  given  in  the  printouts  in  the 
Appendix.  The  errors  that  did  oeeur  are  due  to  a  eombination  of  the  inevitable  eomputational 
errors  and  small  errors  in  the  eurve  fits  to  the  input  MTEs.  Note  that  for  this  case  of  a  sine  wave 
target  input  replaeing  a  bar  target  input.  Equation  (2)  has  a  different  denominator,  given  by: 
f[l-(4Ef)2]. 

•  Processed  bar  targets,  with  various  numbers  of  bars,  with  a  diffraetion-limited  opties  MTE;  the 
latter  is  a  known,  eiosed  form  equation  out  to  the  opties  eut-off  frequeney.  The  output  CTEs 
were  aeeurate  to  within  the  aeeuraey  that  eould  be  extraeted  from  published  data  (Bass  and 
VanStryland,  1994). 

•  Verified  that  the  polynomial  eurve  fits  to  the  MTEs  in  step  (a)  were  smooth,  monotonieally 
decreasing  eurves  with  no  'hidden'  dips  or  spikes. 


15 


3.1.3  Results 


The  CTF  results  for  bar  targets  containing  5,  15,  and  50  bars  are  plotted  in  Figures  4  and  5  for  the 
500  ppi  case  and  1000  ppi  case,  respectively.  Note  that  the  CTF  approaches  a  'steady  state'  as 
the  number  of  target  bars  increases,  e.g.,  there  is  very  little  difference  between  the  15  bar  and  50 
bar  CTFs.  At  the  opposite  extreme,  however,  for  a  low  number  of  5  target  bars,  there  is  a 
distinct  rise  in  the  CTF. 

The  most  relevant  spec  CTF,  for  comparison  to  test  data  results  from  a  real  scanner,  would  be  to 
use  the  spec  CTF  derived  from  the  same  number-of-bars  target  as  is  actually  used  in  the  testing. 
This  is  impractical,  however,  because  of  the  wide  variety  of  bar  target  designs,  both  COTS  and 
specially  fabricated,  which,  for  example,  could  have  a  different  number  of  bars  in  each  frequency 
pattern  in  the  same  target.  Because  nearly  all  practical  bar  targets  would  in  fact  have  between  5 
and  50  bars  at  any  given  frequency,  it  is  convenient  to  define  a  single  spec  CTF  as  the  average  of 
the  5,  15,  and  50  bar  CTFs.  Taking  this  one  step  further,  the  final  spec  CTFs  at  500  ppi  and 
1000  ppi  are  defined  as  the  curve  fits  to  the  average  bar  CTFs,  which  are  shown  in  Figure  6  (see 
also  Table  1).  Note  that,  as  fully  expected  by  the  underlying  theory,  the  spec  CTFs  are 
everywhere  higher  than  their  respective  spec  MTFs. 

One  processing  anomaly  did  occur.  The  CTF  values  calculated  in  the  neighborhood  of  2/3  of  the 
Nyquist  frequency,  i.e.,  close  to  6.56  cy/mm  for  500  ppi  and  13.12  cy/mm  for  1000  ppi,  showed 
an  unnatural  dip  in  modulation;  the  calculated  values  are  only  about  92%  of  what  would  be 
expected  from  a  smooth  curve.  Some  effort  was  expended  in  trying  to  track  down  the  cause  of 
this  but  it  remains  an  unknown.  The  values  near  2/3  Nyquist  are  reported  in  the  Appendix  but 
were  not  used  in  generating  the  spec  CTFs  because  they  are  believed  to  be  anomalous  data 
points,  i.e.,  6.5  cy/mm  for  500  ppi;  12.5  and  13.0  cy/mm  for  1000  ppi  computed  modulation 
values  were  not  used. 
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Modulation 


Spatial  Frequency  (cy/mm) 


Figure  4.  Computed  500  ppi  CTFs  for  5,  15,  50-Bar  Targets  and  Spee  MTF 


Spatial  Frequency  (cy/mm) 


Figure  5.  Computed  1000  ppi  CTFs  for  5,  15,  50-Bar  Targets  and  Spee  MTF 


Spatial  Frequency  (cy/mm) 


Figure  6.  Final  Spec  CTFs  Compared  to  Original  Spec  MTFs 
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3.2  Alternative  Approach  -  Coltman  Formula 

For  analog  imaging  systems,  such  as  a  photographic  film  camera,  it  is  possible  to  measure  the  sine 
wave  MTF  up  to  the  frequency  at  which  the  MTF  goes  to  zero,  which  is  called  the  'cut-off 
spatial  frequency  in  the  optical  field,  equivalent  to  the  'band-limif  frequency  in  electrical 
engineering  parlance.  For  this  case,  a  simple  series  expansion  formula  derived  nearly  half  a 
century  ago  (Coltman,  1954)  will  quite  accurately  convert  the  square  wave  CTF  to  its  equivalent 
sine  wave  MTF,  or  vice  versa. 


Given  the  CTF,  the  Coltman  formula  to  determine  the  MTF,  is 

M(f)  =  - 
4 


C(3f)  C(5f)  C(7f)  C(llf)  C(13f)  C(15f)  C(17f)  C(19f) 

C(f)+ - - - + - + - - - - - - - + - . 


1 1 
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m 


(10) 


and  given  the  MTF,  the  Coltman  formula  to  determine  the  CTF,  is 


C(f) 


4  M(3f)  M(5f)  M(7f)  M(9f)  M(llf)  M(13f)  M(15f)  M(17f) 

—  M(f)  - - + - - - + - - - + - - - + - 

tiL  3  5  7  9  11  13  15  17 


M(19f) 
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(11) 


where,  M(f)  =  sine  wave  MTF 
C(f)  =  bar  target  CTF 
f  =  spatial  frequency 

Note  that  Equation  (10)  is  an  irregular  sequence  in  terms  of  sign  and  term  constants,  whereas 
Equation  (1 1)  is  a  strictly  alternating,  sequential  series.  The  practical  application  of  the  formula 
can  be  illustrated  by  the  following  simple  example.  Suppose  a  CTF  has  been  measured  with  a  bar 
target  out  to  an  imaging  system's  cut-off  frequency  of  30  cy/mm,  then  to  compute  the  sine  wave 
MTF  at  a  single  frequency,  f  =  2  cy/mm,  requires  use  of  the  CTF  modulation  values  at  f  =  2,  6, 

10,  14,  22,  26,  and  30  cy/mm,  i.e., 

.  C(3»2)  C(5»2)  ,  C(7»2)  ,  C(ll»2)  C(13*2)  C(15*2)l 

M(2)--|_C(2)+^  _  +  (12) 


The  fact  that  the  frequency  of  evaluation  increases  with  each  succeeding  term  in  the  expansions 
given  in  Equations  10  and  1 1  indicates  that  for  evaluation  at  higher  and  higher  frequencies,  fewer 
and  fewer  terms  are  used,  because  no  modulation  values  are  available  above  the  cut-off  frequency. 
In  fact,  for  evaluation  frequencies  greater  than  1/3  of  the  cut-off  frequency,  only  the  first  term 
exists,  resulting  in  the  simple  relationships. 


M(f)  =  ^C(f) 


and  C(f)=-M(f)  ,forf> 

7C 


cut  -  off  frequency 
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There  are  several  possible  approaches  for  applying  the  Coltman  formula  to  derive  the  CTF  from 
the  MTF  of  a  sampled  imaging  system,  these  are: 

(1)  Simply  multiply  the  MTF  by  the  4/71  leading  constant  in  Equation  (11);  this  is  a  crude 
application,  but  is  sometimes  seen  in  the  literature. 

(2)  Apply  the  full  Equation  (1 1)  to  the  MTF,  using  as  many  terms  as  is  possible,  given  that  all 
terms  higher  than  the  Nyquist  frequency  are  zero-valued. 

(3)  Extrapolate  the  MTE  curve  beyond  Nyquist  frequency,  until  it  smoothly  reaches  zero 
modulation;  then  apply  the  full  Equation  1 1  with  as  many  terms  as  is  possible. 

Eigure  7  gives  the  CTE  results  of  these  three  approaches  when  using  the  spec  MTE  at  500  ppi  as 
the  input  to  the  Coltman  formula  (Equation  (11)),  and  compares  the  resulting  Coltman  CTEs  to 
the  spec  CTF  derived  via  the  technique  described  in  Section  3.1.  Clearly,  the  4/71  approach  and 
spec-MTF-truncated-at-Nyquist  approach  lead  to  large  errors  at  a  significant  number  of 
frequencies.  On  the  other  hand,  the  spec-MTF-extrapolated-beyond-Nyquist  approach  gives 
results  very  close  to  the  spec  CTF  derived  in  Section  3.1. 

In  fact,  the  integration  kernel  in  Equation  (8)  has  similarities  to  the  Coltman  formula  in  Equation 
1 1,  if  one  expands  Equation  (8)  into  a  series.  To  illustrate  this  similarity,  Eigure  8  shows  the 
CTF  of  a  diffraction-limited  lens  (difflim)  having  a  cutoff  frequency  of  10.0  cy/mm,  as  computed 
by  the  Coltman  formula  and  by  the  method  in  Section  3.1.  For  the  latter,  as  a  further  test  of  the 
computer  program,  the  CTF  for  a  50-bar  target  was  computed  two  ways,  which  should  give 
equivalent  results:  (1)  point  sampling  was  assumed,  with  a  pixel  width  =  0.0001  mm,  pixel  pitch 
=  0.050  mm  (Nyquist=10.0  cy/mm)  and  using  difflimMTF  as  input,  and  (2)  a  finite  pixel 
sampling  width  =  0.050  mm  was  assumed  (Nyquist=10.0  cy/mm,  pixel  pitch=0.050  mm),  using 
[difflimMTF  /  sinc(.057lf)]  as  input.  As  is  seen  in  Figure  8,  all  three  computed  CTF  curves 
overlap,  and  they  also  agree  with  published  data  for  the  CTF  of  a  diffraction-limited  lens  (Bass 
and  VanStry  land,  1994). 

In  summary,  the  technique  described  in  Section  3.1,  which  admittedly  does  involve  more  complex 
computations  than  the  Coltman  approach,  is  believed  to  be  the  more  accurate  approach  for  the 
application  at  hand.  The  Section  3.1  approach  can  account  for  a  small  number  of  target  bars, 
whereas  the  Coltman  formula  inherently  assumes  an  infinite  number  of  bars.  Also,  the  section 
3.1  approach  is  more  conducive  to  modeling  the  discrete-area  sampling  of  a  sampled  imaging 
system,  including  pixel  widths  that  may  be  less  than  pixel  pitch,  whereas  the  Coltman  formula  is 
fundamentally  for  an  analog,  continuous  imaging  system. 
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Spatial  Frequency  (cy/mm) 

Figure  7.  Variations  of  Coltman  Formula  Application 


21 


0123456789  10 

Spatial  Frequency  (cy/mm) 


cy/mm 

difflimMTF 

(input) 

Coltman  CTF 

sec  3.1  CTF 
(point  sampling) 

sec  3.1  CTF 
(area  sampling) 

1 

0.872889 

0.917271 

0.91875183 

0.91875183 

2 

0.74706 

0.830332 

0.8348011 

0.8348012 

3 

0.623838 

0.778428 

0.7798474 

0.77984791 

4 

0.504632 

0.642517 

0.64397621 

0.64397734 

5 

0.391002 

0.497839 

0.49920773 

0.4992092 

6 

0.284757 

0.362564 

0.3637381 

0.363738 

7 

0.18812 

0.239522 

0.239309 

0.2404192 

8 

0.104088 

0.132529 

0.133088 

0.13308893 

9 

0.037386 

0.04760133 

0.04783677 

0.04783712 

10 

0.000000 

0.000000 

0.00015346 

0.00015338 

Figure  8.  Lens  CTF  Calculated  via  Coltman  and  Section  3.1  Method 
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3.3  CTF  Normalization 


In  actual  computation  of  a  scanner's  CTF  from  a  scanned  bar  target,  the  maximum  image 
modulation  is  first  determined  from  the  highest  peak  and  lowest  valley  combination  found  within 
one  period.  It  would  often  lead  to  a  more  accurate  result  if  the  peak  and  valley  image  gray  levels 
were  first  converted  to  their  equivalent  values  in  target  space  before  computing  image  modulation, 
because  this  normalizes  out  gray  level  measurement  differences  between  the  target  and  scanner- 
produced  image.  Bar  targets,  however,  are  rarely  fabricated  with  the  surrounding  uniform  gray 
patches  (see  Figure  1)  necessary  to  construct  the  conversion  curve.  This  leaves  the  choice  of 
either  denoting  the  curve  generated  from  the  image  modulations  as  the  scanner  CTF,  or 
normalizing  this  curve  by  a  very  low  frequency  modulation  value  and  denoting  the  result  as  the 
scanner  CTF.  On  the  upside,  the  latter  approach  has  the  advantage  of  normalizing-out  the  effect 
of  a  bar  target  whose  effective  modulation,  as  'seen'  by  the  particular  scanner,  may  be  less  than 
1.0.  On  the  downside,  the  latter  approach  will  normalize  out  the  degrading  effect  of  any  existing 
light  flare  or  veiling  glare  that  is  a  part  of  the  scanner  performance/design.  Analysis  of  test  data 
from  several  IQS  certified  scanners  that  used  bar  targets,  indicate  that  CTF  normalization  by  the 
modulation  at  a  very  low  frequency  (less  than  5%  of  Nyquist)  is  usually  possible  and 
straightforward,  and  results  in  a  good  measure  of  the  CTF.  Also,  it  should  be  noted  that  if  the 
target  bar  modulations  decrease  with  increasing  frequency,  and  those  target  modulations  can  be 
measured,  then  this  would  be  an  additional  normalization  factor  to  apply  to  the  imaging  system's 
measured  CTF. 
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Appendix 

Computed  CTFs 

Computer  Print-Outs  of  N-Bar  CTFs  Used  to  Generate  Proposed  Spec  CTFs  in  Table  1 
Nomenclature: 

bestphaseCTF  -  used  for  final  spec  CTF;  centers  pixel  on  bar  image  peak  and  valley  for 
maximum  modulation  (see  Section  3.1,  step  (e)). 

PeakCTF  -  maximum  modulation  from  pixels  centered  at  integer  multiples  of  pixel  pitch. 

T(w)  =  equation  4,  5,  6,  or  7,  which  are  curve  fits  to:  (FBI  IQS  spec  MTF)  divided  by  sinc(y). 
5 -Bar  Target  -  implies  there  are  5  target  bars  at  each  target  frequency  (cy/mm)  evaluated. 
Target  frequency  sequence  type  -  defines  the  user-selected  bar  frequency  values. 


5-BAR  TARGET  @  500  PPI 


T (w)  case  (first  term):  4 . 8597D-07*f **7 

tgt  frequency  sequence  type:  1  tgt  form:  BAR 
#tgt  bars  =  5  tgt  mod  =  1.00000000000000 

integral  inc.in  computing  each  BarSinImg  value  =  1.0000 

#  pts  in  integral  to  compute  each  BarSinImg  value  =  4000 

max  frEquation  of  input  MTF  =  40.0001  cy/mm 

pixelwidth  (mm)  =  5.080000E-02  pixel  pitch  (mm) 

#  points  across  pixel  for  pixel  integration  =  401 

pixel  integration  increment  =  1.270000E-04 


OOE-04 

01 


5 


cy/mm 


080000E-02 


cy/mm 

PeakCTF 

bestphaseCTF 

0  . 

.50000000 

0 

.  97679330 

0  . 

.  97679330 

1 . 

.00000000 

0 

.  94837350 

0  . 

.  94837720 

1 . 

.50000000 

0 

.  91613640 

0  . 

.91629111 

2  . 

.00000000 

0 

.  87563520 

0  . 

.  87564660 

2  . 

.50000000 

0 

.  83397450 

0  . 

.  83399051 

3. 

.00000000 

0 

.79415291 

0  . 

.79486203 

3. 

.50000000 

0 

.75908440 

0  . 

.75976770 

4  . 

.00000000 

0 

.71254193 

0  . 

.  72042240 

4  . 

.50000000 

0 

.  67889170 

0  . 

.  68077892 

5. 

.00000000 

0 

.  64079840 

0  . 

.  64088684 

5. 

.50000000 

0 

.59598500 

0  . 

.  60127440 

6 . 

.00000000 

0 

.  54082632 

0  . 

.56083160 

6. 

.50000000 

0 

.46290740 

0  . 

.48141790 

7  . 

.00000000 

0 

.44677844 

0  . 

.48604760 

7  . 

.50000000 

0 

.  42327590 

0  . 

.44946280 

8  . 

.00000000 

0 

.39748110 

0  . 

.41310470 

8  . 

.50000000 

0 

.36846140 

0  . 

.  37629732 

9. 

.00000000 

0 

.33657401 

0  . 

.  33944490 

9. 

.50000000 

0 

.30310210 

0  . 

.30353140 

10  . 

.00000000 

0 

.26907180 

0  . 

.26907180 
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15-BAR  TARGET  @  500  PPI 


T (w)  case  (first  term):  4 . 8597D-07*f **7 

tgt  frequency  sequence  type:  1  tgt  form:  BAR 

#tgt  bars  =  15  tgt  mod  =  1.00000000000000 

integral  inc.in  computing  each  BarSinImg  value  =  l.OOOOOOE-04  cy/mm 

#  pts  in  integral  to  compute  each  BarSinImg  value  =  400001 

max  frEquation  of  input  MTF  =  40.0001  cy/mm 

pixelwidth  (mm)  =  5.080000E-02  pixel  pitch  (mm)  =  5.080000E-02 

#  points  across  pixel  for  pixel  integration  =  401 

pixel  integration  increment  =  1.270000E-04 


cy/mm 
0.50000000 
1 .00000000 

1.50000000 
2 .00000000 

2.50000000 
3.00000000 

3.50000000 
4 .00000000 

4.50000000 
5.00000000 

5.50000000 
6.00000000 

6.50000000 
7 .00000000 

7.50000000 
8 .00000000 

8.50000000 
9.00000000 

9.50000000 
10.00000000 


PeakCTF 
0 . 97703284 
0 . 94693982 
0 . 91353881 
0 . 87199562 
0 . 82957100 
0.78903001 
0.75308910 
0.70610490 
0 . 67195153 
0 . 63345640 
0.58820600 
0.53290070 
0.45529990 
0.43885880 
0 .41498780 
0.38921110 
0 .36079174 
0 . 32971152 
0.29542250 
0.25924620 


bestphaseCTF 
0 . 97703284 
0 . 94694380 
0 . 91369820 
0 . 87200784 
0 . 82958820 
0 .78974200 
0.75379100 
0 .71400433 
0 . 67384120 
0 . 63355034 
0.59353670 
0.55290470 
0 .47380100 
0 .47769320 
0 .44084750 
0.40460040 
0 .36844462 
0 . 33244830 
0.29579880 
0.25924620 
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50-BAR  TARGET  @  500  PPI 

T (w)  case  (first  term):  4 . 8597D-07*f **7 

tgt  frequency  sequence  type:  1  tgt  form:  BAR 

#tgt  bars  =  50  tgt  mod  =  1.00000000000000 

integral  inc.in  computing  each  BarSinImg  value  =  l.OOOOOOE-04  cy/mm 

#  pts  in  integral  to  compute  each  BarSinImg  value  =  400001 

max  frEquation  of  input  MTF  =  40.0001  cy/mm 

pixelwidth  (mm)  =  5.080000E-02  pixel  pitch  (mm)  =  5.080000E-02 

#  points  across  pixel  for  pixel  integration  =  401 

pixel  integration  increment  =  1.270000E-04 


cy/mm 
0.50000000 
1 .00000000 

1.50000000 
2 .00000000 

2.50000000 
3.00000000 

3.50000000 
4 .00000000 

4.50000000 
5.00000000 

5.50000000 
6.00000000 

6.50000000 
7 .00000000 

7.50000000 
8 .00000000 

8.50000000 
9.00000000 

9.50000000 
10.00000000 


PeakCTF 
0 . 98126300 
0 . 94846254 
0 . 91395604 
0.87170910 
0 . 82882070 
0 .78790470 
0 .75167270 
0 .70440262 
0 . 67005930 
0 . 63141750 
0.58603680 
0.53066414 
0 .46432414 
0 .43664741 
0 .41273100 
0.38694770 
0.35864913 
0 . 32747760 
0.29352520 
0.25693010 


bestphaseCTF 
0 . 98126300 
0.94846640 
0 . 91411490 
0 . 87172140 
0 . 82883781 
0.78861860 
0 .75237691 
0 .71232080 
0 . 67195220 
0 . 63151110 
0.59137060 
0.55151730 
0.49905410 
0.47538253 
0 .43851410 
0.40228200 
0.36626520 
0.33020030 
0.29389880 
0.25693010 


50  CYCLE  SINE  WAVE  TARGET  @  500  PPI 

T (w)  case  (first  term):  4 . 8597D-07*f **7 

tgt  frequency  sequence  type:  2  tgt  form:  SINE 

#tgt  bars  =  50  tgt  mod  =  1.00000000000000 

integral  inc.in  computing  each  BarSinImg  value  =  l.OOOOOOE-04  cy/mm 

#  pts  in  integral  to  compute  each  BarSinImg  value  =  400001 

max  frEquation  of  input  MTF  =  40.0001  cy/mm 

pixelwidth  (mm)  =  5.080000E-02  pixel  pitch  (mm)  =  5.080000E-02 

#  points  across  pixel  for  pixel  integration  =  401 

pixel  integration  increment  =  1.270000E-04 


cy/mm 
1 .00000000 
2 .00000000 
3.00000000 
4 .00000000 
5.00000000 
6.00000000 
7 .00000000 
8 .00000000 
9.00000000 
10.00000000 


PeakMTF 
0 . 90816210 
0.79893451 
0 .69179350 
0 .58297440 
0.51372690 
0 .41670331 
0.33931073 
0.30096500 
0.25462052 
0.20003460 


bestphaseMTF 
0 . 90821391 
0.79903161 
0 .69499970 
0.59965870 
0.51386642 
0 .43773480 
0 . 37199431 
0 .31289473 
0.25667250 
0.20003460 
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5-BAR  TARGET  @  1000  PPI 


T (w)  case  (first  term) :  -2 . 99661568D-8*F**5  tgt  frequency  sequence  type 

1  tgt  form:  BAR 

#tgt  bars  =  5  tgt  mod  =  1.00000000000000 

integral  inc.in  computing  each  BarSinImg  value  =  l.OOOOOOE-04  cy/mm 

#  pts  in  integral  to  compute  each  BarSinImg  value  =  800001 

max  frEquation  of  input  MTF  =  80.0001  cy/mm 

pixelwidth  (mm)  =  2.540000E-02  pixel  pitch  (mm)  =  2.540000E-02 

#  points  across  pixel  for  pixel  integration  =  401 

pixel  integration  increment  =  6.350000E-05 


cy/mm 
0.50000000 
1 .00000000 

1.50000000 
2 .00000000 

2.50000000 
3.00000000 

3.50000000 
4 .00000000 

4.50000000 
5.00000000 

5.50000000 
6.00000000 

6.50000000 
7 .00000000 

7.50000000 
8 .00000000 

8.50000000 
9.00000000 

9.50000000 
10.00000000 

10.50000000 
11.00000000 

11.50000000 
12 .00000000 

12.50000000 
13.00000000 

13.50000000 
14 .00000000 

14.50000000 
15.00000000 

15.50000000 
16.00000000 

16.50000000 
17 .00000000 

17.50000000 
18.00000000 

18.50000000 
19.00000000 

19.50000000 
20.00000000 


PeakCTF 
0 . 97750970 
0 . 95400960 
0 . 93073450 
0 . 90736390 
0 .88434980 
0 .86046230 
0 . 83804310 
0 .81657850 
0.79321503 
0.77085030 
0.74581590 
0.72356490 
0 .70161320 
0 . 67917484 
0 . 65614241 
0 . 63364110 
0 . 61787920 
0 . 60188090 
0.58564031 
0.56913950 
0.55219800 
0.53418870 
0.51410690 
0.49102530 
0 .46431261 
0 . 42704750 
0 . 42130103 
0 .41307610 
0 .40473670 
0.39624590 
0 .38754743 
0 . 37857404 
0.36928410 
0.35964372 
0.34961660 
0.33916560 
0.32826894 
0.31692810 
0.30516240 
0.29239472 


bestphaseCTF 
0 . 97750991 
0 . 95401114 
0 . 93073570 
0 . 90736824 
0 .88435852 
0 .86063790 
0 . 83830273 
0 .81659954 
0.79402680 
0.77088140 
0 .74780273 
0 . 72466410 
0.70162590 
0 . 67984990 
0 . 65970870 
0 . 64026212 
0 . 62146550 
0 . 60327830 
0.58586690 
0.56919062 
0.55318641 
0.53781092 
0 . 52297250 
0 .50762522 
0 .48273330 
0 .44431722 
0 .46553412 
0 .45137882 
0 .43725323 
0 . 42317184 
0.40915940 
0.39526300 
0.38159910 
0.36824340 
0.35517620 
0 . 34236764 
0 . 32978913 
0 .31740531 
0.30518963 
0.29239472 
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15-BAR  TARGET  @  1000  PPI 


T (w)  case  (first  term):  -2 . 99661568D-8*F**5 

tgt  frequency  sequence  type:  1  tgt  form:  BAR 

#tgt  bars  =  15  tgt  mod  =  1.00000000000000 

integral  inc.in  computing  each  BarSinImg  value  =  l.OOOOOOE-04  cy/mm 

#  pts  in  integral  to  compute  each  BarSinImg  value  =  800001 

max  frEquation  of  input  MTF  =  80.0001  cy/mm 

pixelwidth  (mm)  =  2.540000E-02  pixel  pitch  (mm)  =  2.540000E-02 

#  points  across  pixel  for  pixel  integration  =  401 

pixel  integration  increment  =  6.350000E-05 


cy/mm 
0.50000000 
1 .00000000 

1.50000000 
2 .00000000 

2.50000000 
3.00000000 

3.50000000 
4 .00000000 

4.50000000 
5.00000000 

5.50000000 
6.00000000 

6.50000000 
7 .00000000 

7.50000000 
8 .00000000 

8.50000000 
9.00000000 

9.50000000 
10.00000000 

10.50000000 
11.00000000 

11.50000000 
12 .00000000 

12.50000000 
13.00000000 

13.50000000 
14 .00000000 

14.50000000 
15.00000000 

15.50000000 
16.00000000 

16.50000000 
17 .00000000 

17.50000000 
18.00000000 

18.50000000 
19.00000000 

19.50000000 
20.00000000 


PeakCTF 
0 . 97774374 
0 . 95255273 
0 . 92808120 
0 . 90367323 
0 . 87971723 
0 . 85494834 
0 . 83171212 
0 .80951330 
0 .78541320 
0 .76240640 
0 .73672430 
0 .71396541 
0 .69158240 
0 .66863930 
0 . 64505141 
0 . 62224274 
0 . 60611402 
0.58975911 
0.57315313 
0.55632940 
0.53909710 
0.52080154 
0.50047051 
0 .47721084 
0.45050784 
0 .41387933 
0.40785700 
0 .39947770 
0.39101650 
0 .38243311 
0.37366540 
0.36456620 
0 . 35512140 
0.34536024 
0.33520520 
0.32462360 
0.31373930 
0.30268010 
0.29139900 
0.27907571 


bestphaseCTF 
0 . 97774392 
0 . 95255450 
0 . 92808234 
0 . 90367811 
0 . 87972690 
0 . 85513270 
0 . 83198370 
0 .80953620 
0.78623360 
0 .76244032 
0.73876351 
0.71506673 
0 .69159320 
0 .66933690 
0 . 64868990 
0 . 62884420 
0 . 60967022 
0.59112620 
0.57336622 
0.55638740 
0 . 54011370 
0 . 52446820 
0 .50937122 
0.49375233 
0.46888460 
0.43104530 
0 .45118070 
0.43693290 
0 . 42274883 
0.40864550 
0 .39464414 
0.38070890 
0.36699384 
0.35361140 
0.34050872 
0.32765520 
0.31515643 
0.30310881 
0.29141830 
0.27907571 
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50-BAR  TARGET  @  1000  PPI 


T (w)  case  (first  term):  -2 . 99661568D-8*F**5 

tgt  frequency  sequence  type:  1  tgt  form:  BAR 

#tgt  bars  =  50  tgt  mod  =  1.00000000000000 

integral  inc.in  computing  each  BarSinImg  value  =  l.OOOOOOE-04  cy/mm 

#  pts  in  integral  to  compute  each  BarSinImg  value  =  800001 

max  frEquation  of  input  MTF  =  80.0001  cy/mm 

pixelwidth  (mm)  =  2.540000E-02  pixel  pitch  (mm)  =  2.540000E-02 

#  points  across  pixel  for  pixel  integration  =  401 

pixel  integration  increment  =  6.350000E-05 


cy/mm 
0.50000000 
1 .00000000 

1.50000000 
2 .00000000 

2.50000000 
3.00000000 

3.50000000 
4 .00000000 

4.50000000 
5.00000000 

5.50000000 
6.00000000 

6.50000000 
7 .00000000 

7.50000000 
8 .00000000 

8.50000000 
9.00000000 

9.50000000 
10.00000000 

10.50000000 
11.00000000 

11.50000000 
12 .00000000 

12.50000000 
13.00000000 

13.50000000 
14 .00000000 

14.50000000 
15.00000000 

15.50000000 
16.00000000 

16.50000000 
17 .00000000 

17.50000000 
18.00000000 

18.50000000 
19.00000000 

19.50000000 
20.00000000 


PeakCTF 
0 . 98197472 
0 . 95407820 
0 . 92849890 
0 . 90339250 
0 . 87891352 
0 . 85371650 
0 . 83011150 
0 .80759763 
0 .78319823 
0.75994170 
0 .73401890 
0.71106112 
0 . 68848431 
0 .66537950 
0.64164102 
0 . 61869070 
0 . 60243880 
0.58596820 
0.56925340 
0.55232830 
0.53500550 
0.51663661 
0 .49624370 
0 .47296091 
0 .44643214 
0 .41675290 
0.40366550 
0.39522720 
0.38673520 
0.37812500 
0 .36931714 
0.36022050 
0.35078150 
0 . 34098420 
0.33082351 
0.32031953 
0.30950270 
0.29838440 
0.28701233 
0.27539940 


bestphaseCTF 
0 . 98197490 
0 . 95407990 
0 . 92850010 
0 . 90339740 
0 . 87892330 
0 . 85390220 
0 . 83038514 
0 .80762080 
0.78402360 
0.75997580 
0.73606914 
0.71216700 
0 . 68849492 
0 .66607950 
0 . 64528860 
0 . 62529772 
0 . 60599600 
0.58733410 
0.56946610 
0.55238664 
0.53602220 
0.52029640 
0.50512051 
0.48998230 
0.47165050 
0 .44792620 
0.44673090 
0 .43244910 
0 .41825400 
0 .40414774 
0.39013590 
0.37622830 
0 .36254671 
0.34915664 
0.33607363 
0 . 32331700 
0.31090320 
0.29880714 
0.28703063 
0.27539940 
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50  CYCLE  SINE  WAVE  TARGET  @  1000  PPI 


T (w)  case  (first  term):  -2 . 99661568D-8*F**5 

Target  frequency  sequence  type:  1  tgt  form:  SINE 

#tgt  bars  =  50  tgt  mod  =  1.00000000000000 

integral  inc.in  computing  each  BarSinImg  value  =  l.OOOOOOE-04  cy/mm 

#  pts  in  integral  to  compute  each  BarSinImg  value  =  800001 

max  frEquation  of  input  MTF  =  80.0001  cy/mm 

pixelwidth  (mm)  =  2.540000E-02  pixel  pitch  (mm)  =  2.540000E-02 

#  points  across  pixel  for  pixel  integration  =  401 

pixel  integration  increment  =  6.350000E-05 


cy/mm 
0.50000000 
1 .00000000 

1.50000000 
2 .00000000 

2.50000000 
3.00000000 

3.50000000 
4 .00000000 

4.50000000 
5.00000000 

5.50000000 
6.00000000 

6.50000000 
7 .00000000 

7.50000000 
8 .00000000 

8.50000000 
9.00000000 

9.50000000 
10.00000000 

10.50000000 
11.00000000 

11.50000000 
12 .00000000 

12.50000000 
13.00000000 

13.50000000 
14 .00000000 

14.50000000 
15.00000000 

15.50000000 
16.00000000 

16.50000000 
17 .00000000 

17.50000000 
18.00000000 

18.50000000 
19.00000000 

19.50000000 
20.00000000 


PeakMTF 
0 .96847522 
0 . 92852234 
0 .89223100 
0 . 85773890 
0 . 82478433 
0 .79162722 
0.76113950 
0.73365200 
0.70256330 
0 . 67867183 
0 . 64637392 
0 . 62459892 
0 . 60400611 
0.57883220 
0.54896580 
0.52103820 
0.50848102 
0 .49408240 
0 .47791761 
0.46006971 
0 .44062793 
0.41969910 
0.39740580 
0.37388831 
0 . 34930911 
0.32384523 
0 .31353271 
0.30688430 
0.29993173 
0.29270310 
0.28523200 
0.27754440 
0.26967042 
0.26163520 
0.25346922 
0.24520080 
0.23686010 
0.22843930 
0.21995791 
0.21141960 


bestphaseMTF 

0 .96848170 

0 . 92856341 

0 .89225083 

0 . 85781264 

0 . 82487094 

0 .79327172 

0.76292520 

0.73376584 

0.70573920 

0 . 67879670 

0 . 65289402 

0 . 62798870 

0 . 60404092 

0.58101210 

0.55886660 

0.53756830 

0.51708430 

0 .49738180 

0 .47843120 

0.46020300 

0.44266820 

0.42580060 

0 .40957440 

0.39358240 

0.37433950 

0.35116913 

0.35061680 

0 . 33725312 

0.32439932 

0.31203890 

0 .30015224 

0.28871881 

0.27772143 

0.26714340 

0.25696951 

0.24718783 

0.23778551 

0.22872020 

0.21997001 

0.21141960 


33 


34 


GLOSSARY 


ANSI 

avg 

CCD 

COTS 

CTF 

cy/mm 

FBI 

lAFIS 

I/O 

IQS 

ISO 

mm 

MTF 

NBS 

ppi 

sinemtf 

tgt 

USAF 


American  National  Standards  Institute 
average 

Charge-Coupled  Device 
Commercial  Off-The-Shelf 

Contrast  Transfer  Function 
cycles  per  millimeter 

Federal  Bureau  of  Investigation 

Integrated  Automated  Fingerprint  Identification  System 

Input  /  Output 

Image  Quality  Specification 

International  Organization  of  Standardization 

millimeter 

Modulation  Transfer  Function 
National  Bureau  of  Standards 
pixels  per  inch 

mitre's  computer  program  to  compute  the  sine  wave  MTF  or  square  wave  CTF 
target 

United  States  Air  Force 
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